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AEROELASTIC STABILITY OF PLATES AND CYLINDERS 

William J. Anderson 
The University of Michigan 

ABSTRACT 

Linear stability criteria are presented for the panel flutter of 
|iin plates and thin-walled cylinders. These structures are exposed to 
fluid flow passing parallel to an outer surface. The expression for 
fluid pressure is simplified in order to emphasize the dynamic proper- 
ties of the systems . The pressures are derived from steady flow rela- 
tions (frequency effects are ignored). An arbitrary spatial phase angle 
is included in the pressure expression. As this phase angle is varied 
in a continuous manner, the fluid flow passes from " subsonic” character 
to supersonic character. The results are useful in classifying several 
types of instability and discussing several pathological cases which 
are usually treated separately. 

The analysis is intended to serve as an aid to understanding the 
mechanism of panel flutter; however, it can be applied directly to 
several problems. It Is accurate for the static divergence and "coupled 
mode" flutter of flat panels in supersonic flow, and also for divergence 
problems wherever experimental measurements can supply the values for 


the necessary aerodynamic parameters. One result is to point out the 
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importance of static instability for flat panels in a transonic viscous 
flow. A second result is to illustrate that the asymmetric divergence 
of cylindrical shells is very sensitive to small changes in the pres- 
sure distribution. 
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Aerodynamic pressure parameter; 

Bh ? /l2(l - v 2 ) 

Airy stress function 
Panel thickness 
Aerodynamic pressure constant, Eq. (5) 

Length of panel 
Mach number 
Axial yave number 
Number of modes 

Axial stress resultant due to initial load 
Circumferential stress resultant due to initial load 
Aerodynamic load 
Integer 

Radius of cylinder 
Time 

Floy velocity 

Panel displacement in transverse direction 

Spatial coordinate, flov direction 

Spatial coordinate 

Kronecker delta 

Amplitude constant 

Angular coordinate 

Eigenvalue 

Eigenvalue 
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Fluid density 
Panel density 

Axial stress in cylindrical shell due to shell motion 
Spatial phase shift 
Frequency, rad/sec 
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Fig. 1. Typical panel flutter problems. 

Fig. 2. Flow over an infinitely long, stationary, two-dimensional 
wavy wall. 

Fig Flow over a two-dimensional flat panel. 

Fig. 4 . Stability boundaries for a flat plate. 

Fig 3* Stability boundaries for a cylinder. 
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1. INTRODUCTION 
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The elastic Instability of thin panels exposed to fluid flow is 
under intensive study at the present time. Typical problems involve 
thin-walled structural elements with one surface exposed to fluid flow 
essentially parallel tc the surface. Figure 1 illustrates the flow 
situation for a flat plate and a cylinder. The usual question of in- 
terest is whether the elastic panels incur divergence (static instability) 
or flutter (dynamic inste.bility) at some value of flow velocity. 

The fluid pressures exerted on oscillating panels are difficult to 
derive in many cases. The role of fluid viscosity, frequency of oscil- 
lation, and panel geometry have complicated panel flutter studies to the 
point where the results are often difficult to understand. 

The present study is based on an intuitive simplification of the 
pressure distribution on the panel. It illustrates the effect of the 
spatial distribution of pressures . The pressures are taken from steady 
flow results and are hence independent of the frequency of oscillation. 

The results are valid only for instabilities occurring at relatively 
low frequencies . 

An approximate solution is required because of the nature of the 
assumptions on the pressures. These assumptions are equivalent to a 
specification of the generalised forces on a discrete system. Galerkin* s 
method is used to pose the eigenvalue problem in matrix form. 
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2 . FLUID PRESSURES 


The pressure expression used in this study is motivated by the so- 
lution for flow over an infinitely long, two-dimenjional stationary wavy 
wall (Fig. 2). For the case of inviscid, isentropic flow, one finds 
that a deflection 


w(x) 


= € sin 


2nx 

£ 



yields a pressure of the form 



where t takes the value 0 for a supersonic flow and jt/2 for subsonic 
flow. The solution is not valid near Mach 1. 

The pressure expression given in Eq. (2) is '’exact” within the 
framework of linearized potential flow for the stationary wall under 
consideration. We will view this expression, however, as an approxima- 
tion which has been provided to describe a given physical situation: 
a panel of finite length with viscous flow effects, real gas effects, 
etc. As an example, for transonic flow, McClure [1] measured pressures 
of the form 


p(x) 
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for a stationary wavy wall. The constants K and ijr are functions of 
Mach number, fluid properties and wavelength. McClure found the ampli- 
tude constant K to be near unity. His measured values of ranged from 
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20 to 45*. We hence see that values of if lying between 0 and 90° 
have physical significance In practical cases. 

Let us consider the pressure expression, Eq. (3) as sufficient for 
our purposes. We will generalize this expression slightly by using sub- 
scripts to show the dependence of the constants K and \|r upon the wave- 
length. For a given deflection of a wall 

w(x,t) = e 1(i)t 2 Bqj sin (4) 

m « 1 ^ L 

one then has a pressure expression of the form 

m 1 1 ■» r K m 008 (% + v) 

(5) 

2L 

Note that each term in Eq. (4) represents a wave with length 

i 

In the following examples, it will be assumed that the constants 
and are known. (This is equivalent to assuming that the generalized 
forces are known for the discrete system.) For example, if slender 
wing ( Acker et) theory were used for supersonic flow over a finite panel, 
Eq. (5) would result with = 1 and « 0 for all m. 

3. FLAT PANEL OF FINITE LENGTH 

Consider the case of a two-dimensional flat panel exposed to 
fluid flow over one surface Fig. 3* The plate is of uniform thickness, 
length L and simply supported a,t both ends. The aerodynamic expression 
of Eq. ( 3 ) will be used to provide fluid pressures above the panel. 
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The fluid belov the panel is at rest and at the same static pressure 
as the upper flow. 

The equation of motion for small deflections of the plate is 


D rr “ K x rf + p s h + »(*>*) ■ 0 

ax* at* 


dx 


( 6 ) 


and the boundary conditions are 


g 2 

v(0,t) « w(L,t) - 2_| (o,t) *. SL" (L,t) « 0 
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The solution is assumed to be of the form 


w(x,t) 
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Galerkin’ s method yields a set of coupled, linear algebraic equations 


of motion 
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E. 

m = 1 


4 N xr 2 
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a m = 0 ( q = 1,2, . . »N) 
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m 2 - q^ if m + q is odd 


if m + q is even 


and is the Kronecker delta. 


This is a linear eigenvalue problem in the eigenvalue It is 


non-Hermitian and hence in general we may have complex eigenvalues 
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The characteristic polynomial is solved for the eigenvalue as a function 
N V L 2 

of A and ~~ — . 

To Interpret the stability of the system, we must remember that 
the frequency of oscillation varies as the square root of the eigenvalue: 

(0 * 


and hence 


v(x,t) 


« e 


ix 1 / 2 * 


The square root must be considered a multivalued function of the com- 
plex variable X, If all eigenvalues X are real and positive, then neutral 
stability results. If X is real and negative, static divergence occurs. 

If X is complex, then flutter occurs. 

Results have been calculated for the stability of a panel with no 
membrane tension (N x * 0). Extensive experience with Galerkin’ s method 
as applied to fourth order differential equations has shown excellent 
convergence when four modes are used. Two-mode, four-mode, and eighth 
mode calculations were used here; the results v?ere found to converge 
adequately. 

The stability boundaries shown in Fig. 4 are from a four-mode 
analysis. For this special case, the amplitude constants and the spatial 
phase shift have been set equal for all modes: 



Kg = K 3 = = K 

t 2 = fj = \ = f 
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As a result, the amplitude constant is easily incorporated into the 
ordinate. The figure hence emphasizes the role played by 

The panel is stable for sufficiently low values of, A, regardless 
of the vralue for t* As A increases, however, the panel becomes unstable 
at some critical value of A, This can be either divergence or flutter, 
depending on the value of 

It is interesting that for ^ » 0 (" supersonic” flow) only flutter 
is possible. (Experimental evidence indicates that this theoretical 
solution is correct for f = 0. ) Also, for \Jr * 90° (” subsonic” flow) 
only divergence is possible. These limiting cases are well known. On 
the other hand, for phase angles ^ between 25° and 90°, one encounters 
divergence first and then flutter. 

The results for small values of say from 0° to 40° are important. 
In transonic flow, for Instance, depends upon boundary layer thickness, 
fluid viscosity, etc. If a given test were carried out for varying 
boundary layer properties, the type of instability might well change 
from a dynamic type to a static type because of this spatial phase shift, 
(it must be remembered that the present analysis cannot predict the single- 
degree-of -freedom type of flutter which often typifies transonic flow. 

On the other hand, this analysis is " exact” for simply supported plates 
which diverge and hence is sufficient to predict static instability.) 

For phase angles near 90°, one finds that increasing dynamic 
pressure causes first a static divergence, followed by dynamic instability 
and finally a static divergence. This might be a confusing factor in 
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some subsonic experimental work; where spatial phase angles might be 
near, but not exactly, 90 c . 


k. ASYMMETRIC FLUTTER OF A CYLINDER OF FINITE LENGTH 


The stabl’’ ' of a finite elastic cylinder Fig. 1 will be investi- 
gated in the same spirit as the flat panel. The shell is of uniform 
thickness and unstiffened. Conventional cylindrical coordinates x, 
r, 0 will be used. Donnell's cylinder equations are adequate to describe 
the deflections of interest here: 


DV 


^.N x ^-^^| + i^E + p 8 h^w + p( X> t) = 0(8) 


ax 2 R 2 d © 2 R dx 2 


8t‘ 


= 0 

R dx 2 


(9) 


The boundary conditions are taken to be the freely-supported case: 


n2 

o w 

V = W = — = a 

bx d x 


(at x = 0, L) 


Again, for a deflection of the form 


w(x,0,t) = 


io)t - . rmrx 

e cos nO sin — - 

L 


the fluid forces will be taken as 


p(x,e,t) = 



a iot 


(cos nS) ijE. cos (22£ + 


If one again applies Galerkin’ s method to the equations of motion 
(8) and (9), one obtains a system of linear algebraic equations 
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where 

^ poha) 2 ^ 

\ = — 35 — 

~ _ pu 2 

A - rf d 

and r)g m is defined as for the plate. 

These equations can be solved for the eigenvalues \ as a function 
of the fluid dynamic pressure ratio % and the phase shift \|/. We will 
consider numerical results for a case corresponding to wind tunnel 
tests carried out by Olson [2]. 

N x = 0 

N e = 0 

R = 8.00 inch 

h = 0.004 inch 

i = 1^.4 inch 

v = 0.35 

n =28 

We will again choose 

^1 = ^2 • tn ~ ^ 

Ki = K 2 = ... K n = K 
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The results for a four mode solution are given in Fig. 5 . Here 
it is seen, as for a. flat plate, that for = 0 only flutter can occur. 
For values of ^ between 60° and 120°, there is an unexpected result. 

The case of static divergence does indeed occur, but at relatively large 
values of ft. In this case, if \Jr is not exactly 90°, then flutter can 
occur at much lower values of A. 

This analysis shows the danger inherent in using an aerodynamic 
theory which predicts that t = 90° exactly. Resulting calculations 
might not reveal a flutter situation which occur at a much lower dyn am ic 
pressure ratio. 

Note that the flutter boundary is very insensitive to changes in 
t from -30° to 60° . This means that the details of the pressure distribu- 
tion on the cylinder are not of much importance in the stability analysis. 
This explains why one of the simplest aerodynamic theories, Ackeret 
theory, can be used with success to predict cylinder flutter which occurs 
at low frequencies [ 3 ]. 


5 . CONCLUSIONS 

The appearance of a spatial phase shift as a free parameter in the 
fluid pressure expression results in some new observations. It il- 
lustrates the change, in a continuous manner, from subsonic (or slender 
body) flow character to supersonic character. The intermediate values 
of the phase angle have physical application to the cases of viscous 
transonic flow over flat plates and supersonic flow over cylindrical 
shells . 
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The analysis is limited to two types of elastic instability: 
coupled mode flutter and divergence. The study cannot predict single 
degree-of-freedom flutter because of the use of steady flow relations 
for the fluid forces . 

Several examples were studied in which the pressure amplitudes 
were identical in all modes and the phase angles \|r m were identical in 
all modes. This case was chosen because of its simplicity. Conclusions 
for the flat plate and the cylinder will be discussed separately. 

The flat plate exhibits both divergence and flutter. For 
one range of the spatial phase angle i|r (-90° to -60°), the plate is 
stable ior all dynamic pressures. For a second range of >|r (-60° to 25°), 
only flutter is possible. Finally, for a third range of \|r (25° to 90°), 
divergence is the critical form of instability, occurring at a much lower 
dynamic pressure than flutter. The stability diagram indicates that 
experiments carried out for certain phase angles might be confusing in 
the sense that different regions of stability and instability could be 
observed in turn as the dynamic pressure is raised. 

Divergence occurs for flat plates at a relatively low value of 
dynamic pressure ratio. As a result, divergence may be a distinct 
problem for the case of viscous transonic flow, where previous pressure 
measurements indicate that the necessary phase shift does occur [1 ] * 

The cylinder example studied was for a particular cylinder geometry 
chosen to match the only successful experiments to date. The cylinder 
exhibits coupled mode flutter over the entire phase angle range of 
physical interest. This flutter boundary is surprisingly insensitive 
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to the value of f. This is fortunate from a practical standpoint. It 
means that coupled mode flutter calculations can be carried out for such 
a shell with less attention paid to the details of the spatial pressure 
distribution. 

The occurrence of divergence for the cylinder is not a simple phe- 
nomenon, In the past, divergence has been predicted for some types of 
cylinders in supersonic flow (where axial wavelengths are long compared 
to circumferential wavelengths). For the cylinder studied here the 
divergence would be of little practical interest. Very small phase 
shifts from f = 90 ° cause flutter to occur at much lower dynamic pres- 
sures than divergence. 

It is not prudent to extend the results of this simple analysis 
too far. On the other hand, it can serve as a qualitative aid to investi- 
gators in panel flutter. There are times when the methods of analysis 
are so cumbersome that one restricts his techniques (or his interest) 
to only divergence or to flutter. It is apparent that one must be 
careful to not overlook one of the possible instabilities. 
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Fig. 2. Flow over an infinitely long, stationary, 
two-dimensional wavy wall. 
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Fig. 4. Stability boundaries for a flat plate. 
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